Linear dynamics of the solar convection 
zone: excitation of waves in unstably 
stratified shear flows 
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Abstract. 

In this paper we report on the nonresonant conversion of convectively unstable lin- 
£SJ ' ear gravity modes into acoustic oscillation modes in shear flows. The convectively 

unstable linear gravity modes can excite acoustic modes with similar wave- numbers. 
The frequencies of the excited oscillations may be qualitatively higher than the tem- 
poral variation scales of the source flow, while the frequency spectra of the generated 
oscillations should be intrinsically correlated to the velocity field of the source flow. 
We anticipate that this nonresonant phenomenon can significantly contribute to the 
' production of sound waves in the solar convection zone. 
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£ ; INTRODUCTION 

The excitation and propagation of waves are important for understanding the dy- 
namics of the sun and stars. It is believed that most of the solar mechanical energy 
is accumulated in the turbulent motions in its convection zone. In the convection 
zone the gravitational stratification drives the convective instability providing the 
dynamical activity of this relatively thin region. The dynamics of the solar con- 
vection is studied to explain many observational features of the Sun. Notably, it 
is thought that the solar acoustic oscillations are excited by the turbulence in the 
convection zone [1-5]. 

Lighthill's ideas of aerodynamic sound generation form the basis of the theoretical 
investigation of the wave excitation in a hydrodynamic medium [6,7]. This theory 
of wave excitation by a free turbulence has been generalized for stratified fluids 
by Stein [8]. From a physical point of view, Lighthill's theory of wave generation 
employs the concept of stochastic excitation of oscillations (waves). In Lighthill's 
theory perturbations are described by an inhomogeneous wave equation, with linear 
terms forming the oscillatory part and the inhomogeneous terms standing for the 
source function. The source terms, which may be classified by their multipole order, 
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are stochastically created by the turbulent perturbations. The amplifying effect of a 
sheared mean flow on the fluctuations of the Reynolds stress (nonlinear source term) 
and thus on the wave production has been noted by Lighthill [7]. However, this 
effect has not received further attention within the context of stochastic excitation. 

Significant advances in the investigation of the dynamics of flows with velocity 
shear have been achieved together with the disclosure of specific features of shear 
flow phenomena [9,10]. Operators arising in the mathematical formalism of the 
canonical modal analysis in the study of the linear dynamics of shear flows are not 
self-adjoint. Consequently eigenmode interference introduces principal complica- 
tions. The nonmodal approach has proved to be an alternative successful route 
for exploring the dynamics of shear flows. This approach employs the study of 
temporal evolution of the spatial Fourier harmonics of perturbations. 

Impressive progress has been made by use of the nonmodal analysis (see e.g., [11- 
16]). This approach has led to the discovery of new channels of energy exchange 
between different modes in shear flows. Resonant phenomena of wave transforma- 
tions have been studied in [17-23]. The nonresonant phenomenon of the conversion 
of vortices into acoustic waves has been described in [24]. The same mechanism is 
found to operate for magnetosonic [25] as well as for plasma Langmuir oscillations 
[21]. 

In this report we introduce a new dynamical source of acoustic waves in unstably 
stratified shear flows. Namely, the linear nonresonant conversion of convective into 
acoustic wave modes in a stratified shear flows. Convectively unstable exponentially 
growing buoyancy perturbations generate acoustic wave oscillations in presence of 
a sheared mean flow. We identify this linear conversion of modes in shear flows as a 
new excitation mechanism of the solar oscillations and waves. It differs in principle 
from the stochastic excitation mechanism and should significantly contribute to the 
process of acoustic wave generation in the solar convection zone. 



PHYSICAL APPROACH 

The equations governing the dynamics of a compressible stratified flow are: 

[«9 t + (VV)]p + p(VV) = 0, (l.a) 

[3t + (VV)]V = -VP/p + g, (1.6) 
[d t + (VV)] P = ( 7 P/p) [dt + (VV)] p. (l.c) 

We consider the hydrodynamic situation where a horizontal shear flow V = 
(Ay, 0,0) occurs in a vertically stratified medium g = (0,0,— g). For simplicity 
we assume that A = const and g = const. This yields the stratified equilibrium 
state: 

P (z)/P (0) = p (z)/p o (0) = exp(-zk H ), (2) 



where ku = 19 /c 2 and c 2 = 7P0/P0 We introduce the linear perturbations in the 
following way: 

V = Vo + V'p (0)/po(z), P = P + P', P = Po + P- (3) 

Here the velocity perturbations are normalized to exclude the exponential height 
dependence due to the vertical stratification of the background flow. We use the 
Cowling approximation [26] and neglect the perturbations of the gravitational ac- 
celeration. Following the standard method of nonmodal analysis (see [27] for a 
rigorous mathematical interpretation) we introduce the spatial Fourier harmonics 
(SFH) of the perturbations with time dependent phases: 

*(r, t) = ip(k(t),t) exp(ik x x + \k y {t)y + ik z z), (4. a) 

ky(t) = ky(0) ~ Ak x t, (4.6) 

where k z = k z + ifc#/2. For compactness of notation we introduce the generalized 
vector of perturbations and their SFHs as follows: ^ = (V',p r , p') and ip = (u,p, p). 
To avoid complex coefficients in the dynamical equations, we construct the normal- 
ized entropy and vertical velocity perturbation SFHs in the following way: 

s = (ic 2 s ~k*Jg-l)(p-c 2 s p)/^-l), (5.o) 

v = (c 2 s k* z + ig)u z , (5.6) 

where k* z = k z — ikjj/2. From Eqs. (1-5) we obtain, by the use of straightforward 
manipulations, the following set of differential equations that govern the SFH of 
the linear perturbations in stratified shear flow: 

p(t) = c 2 s (k x u x + ky(t)u y ) + v, (6. a) 

u x (t) = -Au y - k x p, (6.6) 

Uy{t) = -k y (t)p, (6.c) 

v(t) = (N 2 - c 2 s k 2 z )p - N 2 s, (6.d) 

s(t) = v. (6.e) 

is the square of the frequency of the Brunt- Vaisala: = gkuid — l)/7 and 
k 2 z = \k z \ 2 = k 2 + k 2 H / A. In an unstably stratified flow negative buoyancy (iV| < 0) 
requires that the adiabatic index 7 < 1. Such an effective value may be assigned 
to this parameter under a certain thermodynamic approach [28]. However, in Eqs 
(6.a-e) we retain only 7V| and argue that these equations are more general than the 
underling 7 prescription. 

Further we note that vorticity is conserved in the wave-number space: / = 
k x u y — k y (t)u x — (A/c 2 )(p — s). The spectral energy of the perturbations can be 
defined as follows: 

E = p Q /(2c 2 s ) (E K + E P + E T ) , (7. a) 



EK = c 2 s (ul + ul)+v 2 /(cm-Ni), (7.6) 
E F =p 2 , E T = Nls 2 /{c 2 s k 2 z -Nl). (7.c,d) 

where E^, Ep and E^ correspond to the kinetic, elastic and thermobaric energies 
of the perturbations, respectively. Formally the perturbation energy is conserved 

in the shearless limit: E = Ac 2 u x u y . The instability of the convective eddies 
corresponds to a negative value of the thermobaric energy. 



Linear modes in the shearless limit 

The linear modes may be classified explicitly in the shearless limit (A = 0). In 
this case the full Fourier expansion of the linear perturbations \I/(r, t) oc ^(k, uj) 
yields the dispersion equation: 

cu(u 4 - c 2 s k 2 u 2 + N 2 c 2 s k 2 ± ) = 0, (8) 

where k\ = k 2 + k 2 and k 2 = k\ + k 2 . The solutions of Eq. (8) describe the stability 
and characteristic temporal variation scales of the existing modes: 

u v = 0, (9. a) 

where the subscripts v, s, c define the frequencies of the vortex, acoustic and 
convective modes, respectively. In an unstably stratified flow, i. e., when iV| < 0, 
iu c defines the growth rate of the buoyancy perturbations. 

Obviously the / = constant law demonstrates the existence of the stationary 
{uj = 0) vortex mode in the linear spectum. The conserved vorticity I may be 
considered as the vortex mode measure. The physical eigenfunctions of the acoustic 
<3> s (t) and convective $ c (^) modes may be rigorously defined in this limit: 

$ s (t) ee p(t) + K^f^- (s(t) - k f 2 p{t)^ (lO.a) 

Mt) = s(t) - f 2 p(t) - ^^P(t) (10.6) 

where Q 2 . = c 2 k 2 and Vt 2 = N^k\/k 2 . Hence the equations governing the dynamics 
of the perturbations of the different modes may be decoupled as follows: 

'6 s (t)+tu 2 $ s (t) = 0, (ll.a) 

$ c (t) + cu c 2 $ c (t) = 0. (11.6) 

Starting from this simple situation we study the velocity shear effects on the per- 
turbation modes. 



Effects of a sheared flow 



To study the effects of the velocity shear on the linear modes we introduce the 
small-scale approximation: k 2 ^> k 2 H . This approximation strongly simplifies the 
mathematical formulation and is justified for the following two reasons. Firstly, 
our analysis needs constant vertical gravity, an assumption that may be adopted 
for perturbations with vertical height scales shorter than the stratification scale. 
Secondly, this approximation is necessary for our assumption of a constant linear 
shear of the flow velocity, especially in the turbulent flows. Using Eq. (7) this 
approximation may be represented by the following condition: c 2 s k 2 z ^> iV|. In terms 
of the frequencies it yields (Q 2 — uo 2 ) pa (Q 2 — uj 2 ) pa 0, which strongly simplifies 
the characteristic physical quantities of the perturbation modes: $ s (t) ~ p(t) and 
$ c (t) pa (s(t) — k 2 p(t)/k 2 (t)). To analyze the dynamics of acoustic oscillations in the 
shear flow we rewrite Eqs. (6.a-e) in the form of coupled second order differential 
equations for the variables pit) and y(t): 



pit) + f(t)p(t) + n?(f)p(f) = Ax(*)y(0 + A 2 (f)y(f), 
where the convection variable y(t) is introduced as follows: 



y(t) 



kj_(t) 
k(t) 



s(t) - 



kl 



k 2 (t) 



p(t) 



and 



ili(t) = c 2 s k 2 (t) + 2A 



m = N^ + 2A 2 



k 2 k 2 



k\(t)k\t) 



fit) = 2A 



Xi(t) = -2A 



k 2 (t) ki(t)k\ty 

k x ky (t) 



k 2 (t) ' 

k x ky(t) 



k ± (t)k(t)' 



\ 2 {t) = -2A- 



,k 2 x k ± (t) 



k 2 Jt)k 2 2 



(12.a) 
(12.6) 

(13) 

(U.a) 
(14.6) 
(14.c) 
(U.d) 
(14.e) 



k 3 (t) \ k A ± (t) J ' 

In deriving Eqs. (12.a-b) we have used the following two simplifications. Firstly, we 
have retained only the terms describing the effect of the buoyancy perturbations on 
the acoustic waves, and we have neglected the effect of the acoustic pressure pertur- 
bations on the evolution (exponential amplification) of the buoyancy perturbations 
in the right hand side (rhs) of Eq. (12. b). Secondly, we have neglected the source 
terms in the rhs of the two dynamical equations that describe the shear induced 



coupling between the vortex and acoustic wave modes (in Eq. 12. a) and vortex 
and buoyancy modes (in Eq. 12. b). In fact, the coupling of the vortex and acous- 
tic wave is a process that has been studied to reveal the mean flow shear induced 
nonresonant mode conversion phenomenon in [24]. However, in the present case, 
the source terms of the acoustic waves that are proportional to the vortex mode 
measure, conserved quantity /, are dominated by the source terms, associated with 
the exponentially amplifying convective modes: y(t) and y(t). It should be empha- 
sized that the present approach is justified only for a convectively unstable medium 
with iVg < 0, so that the buoyancy modes undergo exponential amplification in the 
linear regime. 

The dynamics of the acoustic waves in the absence of the buoyancy perturbations 
is described by the homogeneous part of Eq. (12. a). The acoustic wave frequency 
and amplitude variations are described by the parameters £li(t) and f(t) (see [19] 
for a detailed study). The dynamics of the convective mode is described by Eq. 
(12. b). Eq. (16. b) shows the transient stabilization effect of the sheared mean 
flow in an unstably stratified medium. The stabilization occurs at times, when 
\k y (t)/k x \ < 1 and reaches its maximum at t — t*, when k y (t*) = (see Eq. 16. b). 

The terms Xi(t)y(t) and X 2 (t)y(t) in the rhs of Eq. (14. a) describe the coupling 
between the convective and acoustic waves modes. The shear flow origin of these 
source terms is obvious from Eqs. (16.d,e). Hence, Eqs. (14.a,b) describe the 
mean flow shear induced buoyancy - acoustic wave mode conversion in a convec- 
tively unstable medium. Some specific features of this phenomenon are due to its 
linear nature; SFH of the exponentially growing buoyancy perturbations are able 
to generate SFH of the acoustic waves with the same wave-numbers. The ampli- 
tude of the excited wave mode depends on the values of the source terms Xi(t) 
and X2{t). So, convective modes with k x = can not generate acoustic waves at 
all (Ai = A2 = 0). While maximal efficiency of the mode conversion phenomenon 
should occur at k z = 0, or in a realistic physical approximation (see Eq. 12) at 
k 2 z > N^/c 2 s . Naturally, acoustic wave emission from convection should generally 
increase when the mean flow shear parameter A increases. 

We numerically analyze Eqs. (6.a-e) to verify the analytical results obtained 
from the approximate equations (12.a,b). We select the initial perturbations in a 
specific manner, which enables us to excite the convective and acoustic wave modes 
individually at the initial moment of time. It appears that exponentially growing 
buoyancy perturbations instantly excite the acoustic wave mode harmonics at a 
given point in time, when the perturbation wave-number along the flow velocity 
shear is zero: t = t*, k y (t*) = 0. The generation of acoustic waves is clearly 
traced from the pressure variation, as well as the compression of the perturbations. 
Numerical analysis shows that the efficiency of this mode conversion phenomenon 
increases with the flow shear parameter. 



DISCUSSION AND CONCLUSIONS 



We have presented a study of compressible convection in shear flows. In particular 
we have focused on linear small-scale perturbations in unstably stratified flows with 
constant shear of velocity. The linear character of the system enables us to identify 
the perturbation modes and to study their dynamics individually. We find a mode 
conversion that originates from the velocity shear of the flow: exponentially growing 
perturbations of convection are able to excite acoustic waves. This process offers a 
novel approach to the hydrodynamic problem of the acoustic wave generation. 

This wave excitation phenomenon can be important for the acoustic oscillations 
of the sun. Being responsible for the wave generation in high shear regions of 
a stratified turbulent flow, this nonresonant phenomenon can contribute to the 
production of sound in the solar convection zone. Moreover, the process of the 
wave excitation should be triggered by a weak vertical magnetic field. In this case 
we anticipate the production of high frequency compressional MHD waves. The 
latter process will considerably increase the extraction of the mechanical energy of 
the convection by waves. 

Specific to this phenomenon is that perturbations of buoyancy are able to excite 
acoustic waves with similar wave-numbers. This property makes it clearly dis- 
tinct from stochastic excitation, where the generated frequencies are similar to the 
life-times of the source perturbations. In contrast, frequencies of the oscillations 
generated by the mean flow velocity shear induced mode conversion may be quali- 
tatively higher than the temporal variation scales of the perturbations in the source 
flow of a compressible convection. The frequency spectrum of the excited acoustic 
waves should be intrinsically correlated to the velocity field of the turbulent source 
flow. Shear flow induced wave excitation in stratified flows offers a natural expla- 
nation of the fact, that the solar acoustic oscillation are mainly excited in the high 
shear regions of the convection, intergranular dark lanes [29]. It also explains the 
puzzling wave-number dependence of the observed mode energies at fixed frequen- 
cies (see [5] and references therein). A detailed comparison with observational data 
requires a more realistic physical model. The simplicity of our model is used to 
demonstrate the basic features of this excitation phenomenon. 

Finally we note that in the present formalism we have focused on the waves 
with frequencies higher than the characteristic cut-off frequency for the acoustic 
waves in the convection zone. Shear flow initiates the qualitative change of the 
temporal variation scales of perturbations and the excitation of the waves that are 
not trapped in the convective envelope. Hence, this mode conversion presents a new 
significant contribution into the channel of energy transfer from the dynamically 
active interior to the atmosphere of the Sun. 
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